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■ Abstract. Chiral symmetry, fundamental in the physics of graphene, guarantees the 
[ existence of topologically stable doubled Dirac cones and anomalous behaviors of the 

■ zero-energy Landau level in magnetic fields. The crucial role is inherited in the optical 

responses and many-body physics in graphene, which are explained in this paper. 
We also give an overview of multilayer graphene from the viewpoint of the optical 
properties and their relation with the chiral symmetry. 



1. Introduction 

While tlie Dirac cone physics in graphene is forming a new branch in the condensed- 
matter physics, there is an increasing fascination with optical responses in graphene [Tl|2l 
[3]. The present article focuses on how the physics of massless Dirac cone and the optical 
properties of graphene are related. There, we stress the factor that links the two is the 
chiral symmetry. So the purpose of this paper is, first, to discuss the meaning of the 
chiral symmetry in terms of fundamental ingredients in graphene including interaction 
and multi-layer systems as well as optical responses. 

So let us start with the Dirac cone, which comes from the symmetry of graphene. 
The honeycomb lattice structure of this material has two atoms per unit cell, and if 
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we focus on the two 7r-bands arising from this, the system is effectively described by a 
traceless hermite Hamiltonian, H{k) = R{k) ■ a, which is parametrized by R with three 
real components. This implies that the Dirac cone has a co-dimension of three Since 
the momentum in two dimensions are described by two parameters, the Dirac cones do 
not arise in general, while the symmetry of the graphene at K and K' points gives rise 
to the Dirac cones. Here the key ingredient is the chiral symmetry of the honeycomb 
lattice. When the lattice points are classified into two sublattices and the hopping of 
electrons is only allowed between them, the system is chiral-symmetric. Then we can 
show that the degree of freedom in R is reduced from three to two. This also ensures 
the topological stability (against various modifications) of the Dirac cones of the chiral 
symmetric two-dimensional system, since the degeneracy cannot be removed by small 
but finite modifications as far as the chiral symmetry (or its extension) is respected. 

A conspicuous feature is graphene and related models is that the Dirac cones 
appear in pairs (at K and K' points in graphene). This is due to the well-known 
fermion doubling in chiral-symmetric systems, which is a two-dimensional analogue of 
the Nielsen-Ninomiya theorem in four-dimensional lattice gauge theory [5l |6]. In a real 
graphene, the chiral symmetry is not rigorously valid. In the standard tight-binding 
parameters for graphene [7], however, most of the chiral-symmetry-breaking parameters 
are significantly smaller than the chiral symmetric ones in magnitude. 



2. Chiral symmetry 



Let us start with the chiral symmetry of fermions on a lattice, which is directly applicable 
to graphene. Assuming the lattice is bipartite, that is, all of the lattice points can be 
divided into two sublattice sites • and o, the non-interacting Hamiltonian "Hq is expressed 
as block-off diagonal form. 



■Ho{D) = c^Hqc = c\Dco + h.c. 
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where N, and A^o are the numbers of • and o sites, respectively, and D is an A^, x A^o 
matrix. We can introduce the chiral operator F that has 

In. O 

^ O -In. 

where = J^, A^ = A^. + A^, and Tr T = A". - A^,. 

The chiral symmetry {H, F} = implies that, if <fi is an eigenstate of H with 
energy Aj Tipi is another eigenstate with energy — Aj {HoTipi = —THoipi = —XiTipi). 

Without loss of generality, let us assume A", > A'o- Then the secular equation 
becomes 



= detjv(A/jv - Ho 
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= A^--^°det;v.(A^/iv.-i^^^). (5) 

This implies, at least, N, — No of the e^'s are zero. These zero energy states (zero modes) 
are geometical since the number of these exact zero-energy states is determined by the 
geometrical structure of the lattice. 



2.1. Chiral basis for the zero modes 



As for a zero-energy state, its chiral partner is degenerated in energy. This implies that 
there is a gauge freedom in choosing the basis. In this paper, we focus on the zero-energy 
Landau level in graphene. 

Before proceeding, we have to take care a bit about what we really mean by the 
zero-energy Landau level. Rigorously speaking, the n = Landau level has the energy 
exactly equal to zero only for massless Dirac fermions in a continuous model. Conversely, 
if you consider a lattice model such as honeycom tight-binding model, the n — Landau 
level has in general a finite width, and the zero-energy condition is only realized in the 
limit of weak magnetic fields (0 — > where is the flux per hexagon). In a numerical 
calculation for the honeycomb lattice model with a finite 0, the n — Landau level has 
a small but finite width. Still, we want to project the problem onto the n = Landau 
level. For this purpose, we can define the e-zero mode by collecting those states that 
have energies between — e and e, where e(> 0) means the width of the Landau level. 
These naturally include the geometrical zero modes discussed above. 

Now let us take an orthonormalized basis for the e-zero modes which form an M- 
dimensional multiplet as ^Z' = (^i, ■ ■ • , tpM) {i^li^j = hj)- 

{4\Hl\i^.) < (6) 



for ^i. This property is inherited by a unitary-transformed ip^ = {{ipuj)i, ■ ■ ■ , {'^ui)m) — 
ipuj for u G U{M). It imlies the e-zcro modes arc U{M) gauge invariant. 
Then one can construct a normalized complete multiplet as 

V'T = (V',¥',r(^)- (7) 

where </? is a multiplet of negative energy states i^oV'i — \ > ^ and Tip is a 

multiplet of the chiral partners of (p. The normalization and the orthogonality of all the 
states are expressed as 

ijj^ ijj^ip ij^Tip 1 r /m " 
ip^i/j ip^ip p^Vp = Im' (8) 
p^Fip p^Fp p^p Im' 

where M' is the dimension of the negative (positive) enerygy multiplet p (Lt/?). The 
completeness reads 

V'tV't = V'V'^ + <ff^ + I^ff^I" = In- (9) 
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With this multiplet ipT^ Hq is block diagonahzed as 

S 
-A 
A 

£ = i)^Hoip, ^^ = Tr^V^<e^ A = diag (Ai, ■ ■ ■ , Aa/') (^i > e). (11) 
Due to the completeness and orthogonality relations Eq.([8]), we have 

T^l) = i/jT^^r^Tip = ifj^^T^ + ifiip^T^ + Tifiip^T'^ilj = i/jTo, (12) 
To = i^^Tij, (13) 

where the M x M-dimensional hermitian matrix Fq satisfies Fq = Im- Its trace is 
evaluated as 

TrFo = TrF(V^V^t) = Tr F - 2Tr F(^^^) = Tr F = A^. - A^^ (14) 
Then Fq is diagonahzed as 

(15) 

/m+ O 
O -Im. 



M+M- 



M+-M_ = N,- N^, 



(16) 



where ur € U{M) with M+ + M_ = M. Now let us define a chiral multiplet ipr = ^i^r 
that satisfies 



F?/>_ = -?/>_. 4)_ = ■ ■ ,^M--) 

Namely, V^j^'s are eigenstates of F as 











(17) 
(18) 
(19) 

(20) 
(21) 



where ip, is an A^, x matrix and ipo an A^o x matrices. 

3. Electron-electron interaction and chiral condensate in a magnetic field 

Let us move on to discuss the electron-electron interaction in graphene in magnetic 
fields. 

3.1. Two-body interactions 

As for the electron-electron interaction, we take the two-body interaction VijUinj 
between sites (ij), which can include long-range interactions. Note that 



TliThj Cj^CiCjCj C^CjCjCi 



-1 + Hi + Uj + CiCjC^jc\. 



(22) 



Chiral symmetry and its manifestation in optical responses in graphene 



5 



We can express the interaction in an electron-hole symmetric form, at half-filling, as 

2 7 V ^ ~ 2 i 2 



H 



int 



ij) 



rii 



+ c 



ct) 



const. 



(23) 



3.2. Ground state within the projected subspace 

We assume that the electron-electron interaction is sufficiently weak as compared with 
the Landau gap between the n — Landau level and those forn = ±1. When the fiUing 
factor of the n = LL is z/ < 1, the ground state is then given by the configurations of 
the n = LL states, while the filled states below them can be regarded as the "Dirac 
sea" . The one-particle spectrum in graphene is given by the e-zero modes and the chiral 
pairs of non-zero energy states as discussed above. Then the kinetic energy from the 
e-zero modes is negfigible, so that the ground state is given peturbatively (i.e., in the 
interaction energy that is smaller than the Landau gap) as 

M' 



n, 



■,iM<Z{l,-,M} 



<h 



l^<)=n^<]|0)' (24) 



where |D<) is the filled Dirac sea ( di\D^) — for all i's), {ii, • • • , im} a subset of the 
e-zero modes, and M the number of particles. The free-fermion Hamiltonian 1-Lq (i.e., 
the kinetic energy) is written as 

Uq = cVt diag {E, -A, k)'^\.c = SSd - d<^ kd< + d>Ud>, (25) 



(26) 



where d^ — (o?|, • • • , o?m), d< = ((i<J, • • • , (i<jv^/), and (i>^ = {d>\, • • • , d>l^,) are standard 
fermion creation operators for the e-zero mode, the negative energies and positive ones 
as {d£,dl} = 6ie, {d<:i,d^^^,} = du>, {d>£,dl^,} = dw, etc. 
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3.3. Projected interaction and its symmetry 

Now let us take the chiral basis for the e-zero modes to describe the interaction between 
the particles. We assume that the perturbative ground state within the configuration of 
the n = Landau level is determined by the projected interaction 1-Lint ■ The projected 
Hamiltonian is defined as 

'Hint ^\Y. V,^[c\cycfc, + (c ^ ct)] = 1 5] Vi^{dH^UdH^)^{il^d)^{il^d)i + C.C., (27) 

ij ij 

where C.c. is a charge conjugation defined below and q is a projected fermion defined 
as 

^ d 1 \ d 

*T =(V',¥?,r(^) 





-0^ = '0'0tc — Pc, 



(28) 
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with P = ipip^ = being the projection to the e-zero modes. These q's are not 
canonical fermions, since 

{q, c]} = {{P)u'C,, 4(pt) ., .} = (ppt) .. ^ (p2)^. ^ (p)^._ (29) 

Note that the projected Hamiltonian is clearly positive semi-definite. 



3.4- Chiral transformation and charge conjugation 

Let us take a chiral basis to define canonical fermions di for the e-zero modes, defined 
as 

d^ 
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We can define a chiral transformation as 



(30) 
(31) 
(32) 

(33) 



where the generator of the transformation, the chirality is given as 

g = c^Tc = J2 4ci - Yl ^ d^i^lrjjjrd = d^Tod = dld+ - dld^. (34) 



The chiral transformation operates as 



d^ de=UldUe = e'"^ de 



e-'^d+ 



Ci ^ 



(35) 
(36) 



Then it is clear that the projected interaction is invariant under the chiral transformation 
as 



'Hint ^ 6 ^Hint^~^ ^ — "Hint- 

For an infinitesimal transformation this implies a conservation law, 

[G, -Hint] = 0. 

Let us next define a charge conjugation which is anti-unitary as 



O ^ Al^OAc, 



M 



M' 



Ac^K Hide + dl) J](d<, + <,)(d>, + <,), 



(37) 

(38) 

(39) 
(40) 



£=1 
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where K is complex-conjugation with K"^ — 1. Its operation is, for example, 

A^cdtAc = H^'-^dl A^'cdlAc = {-f-'de. (41) 
Then the invariance of the interaction Hamiltonian under the charge conjugation, 

A^C^intAc = 'Hint (42) 

trivially follows. 

3.5. Chiral condensate as doubly-degenerate ground states 

In this section, we restrict ourselves to consider a repulsive interaction that only acts 
between • and o sites. The projected interaction in this case is written as 

Tiint ^\Y1 (4c]oCjoCi. + C.c) (43) 

i'jo 

= ^ E [{dUlUdUl)j{^od-)AMi + (MiiMAdUlWUl^ m 

i.jo 

where ^Vi,j^ > 0. Here we allow any repulsive interaction between the • and o sites. 

Now we can readily see that the two states, 

nint\G±) = 0, IGi) = 4± • • • 4±±l^<), (45) 

which have maximum and minimum chiralities, respectively, are special, since each of 
them has vanishing interaction energy. As for for example, all de^ states are filled, 
{dl-^l)i\G+) = 0. Also none of the d^^ states are occupied, {ipod-)j\G^) = 0. Since 
the interaction Hamiltonian 'Hint is semi-positive definite, these two states are ground 
states of Tiint with a degeneracy of two . 

As for the chiral transformation, we have 

|G±)^Wj|G±) = e^^^±^|G±), (46) 

where U^, is the chiral operation defined above. This property follows, since \G±) are 
eigenstates of the chirality Q as 

g\G±) = ±M±\G±). (47) 

Thus we have shown that the ground states are the chiral condensates, where the 
chirality is macroscopic (M± ~ A^). 

As for the charge conjugation, it operates as 

|G±)^(-)M^---4^^|D>), (48) 

where the \Dy) is the positive Dirac sea up to the sign. 

3.6. Chiral-condensate doublet as a Hall insulator 

The Hall conductance of the doubly-degenerate chiral condensates can be calculated 
with the Niu-Thouless-Wu formula as 

o2 



^-^C, C=-^/Tr2F, F^dA + A\ A = *trf*, (49) 
h Nd 2m J 
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where \1/ = (|G'+), is the chiral-condensed ground states with the degeneracy of 

Nd = 2. 

The doublet, being degenerate, can be mixed, but diagonahzation of the chirahty 
Q acts to fix the gauge. Here we assume a finite energy gap above the ground state 
multiplet. Then the Chern number of the doublet becomes well-defined, and given by 
the sum of each as 

C = C+ + C., (50) 

where C± is the Chern number for \G±). Further, we can decompose the condensate as 

C± = + Cd^, (51) 

= ^ y TrAf±rfV±#±, C'd< = (rf/^<MD<), (52) 

where Cd^ is the Chern number of the filled Dirac sea. Since the chiral operator F is 
unitary, the Chern number of the positive Dirac sea is the same as that of the negative 
Dirac sea as 

Ci)> = ^ y {dDy\dD^) = Co^ = Cd. (53) 
Since the charge conjugation is anti-unitary, it implies 

Cd>+Q_ = -(Cd<+C^J, (54) 

which in turn gives the total Chern number of the doubly degenerate chiral condensate 
as 

C = {Cd^ +C^_) + {Cd^ + C^J = + C^^ + 2Cd = 0. (55) 

This implies that the half-filled state formed by the chiral doublet is a Hall insulator 
with the topological degeneracy of 2. 

4. Multilayer graphene and chiral symmetry 
4.I. Bilayer/ trilyer Hamiltonians 

Let us start with the bilayer graphene in magnetic fields[6]. With j = (ji, j2) denoting 
two-dimensional coordinates, ei^2 the unit translations, and if the total flux per hexagon 
in units of the flux quantum (fQ = h/e (see FigH]), the Hamiltonian is written as 

^d = 7o5^[c^l(j)^o(j) 

j 

+e'2-^^Ml(j)t^o(j - 62) + e-*2Mii+|)4(^- + e^)4(^-) + h.c.] , (57) 

nu = Hd{d ^ M, 31 ^ 31 + 2/6) (58) 

y.ud = [liul{3)do{3) 
j 

+73(^e*2.c,o,+i)^t(^-) + di{3 + ei) + e-^^^^^^'^+t^ciKj + ei - e2))^io(j - 62) + h.c] (59) 
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Figure 1. Lattice structure of the bilayer graphene. Interlayer couphngs 71 and 73 
are denoted by green and yehow bonds, respectively. 



Here "u{d)" denote up (down) layers, 70 is the nearest-neighbor hopping within each 
layer, and 71 the inter-layer hopping perpendicular to the layer. The hopping 73, 
the second largest inter-layer hopping, is along an oblique direction, causes trigonal 
deformation of the energy dispersion (trigonal warping) [8] in zero magnetic field. In a 
finite magnetic field, the hopping acquires Peierls phases appearing in the Hamiltonian 
above, again due to the oblique directions of the hopping. 
As for the trilayer, we have 

^trilaycr = + + + + V,t^^'^^^ (60) 

Ht = 'Hd{d^t,3i-^ ji + A/Q) (61) 
H^f"^ = Undid ^u,u^t, Jl ^ Jl + 2/6) (62) 
H:^/''':Hudid^t,u^u). (63) 

Here we have added "t" (top layer) on top of the "d,u" (down, up) layers, and ABA 
(ABC) stand for the trilayer graphene with the ABA (ABC) stacking. 

Now let us Fourier-transform the fermion operators, first along the 62 direction as 

-^e''"'^''da,Ua,ta{ji,k2). (64) 



Then we have 



= / -^n'^^'^ih), = / ^n'^'^h), (65) 



n'^^'^ih) = nY'ih) + nlTih) + nl^ih), (66) 
n'^'^h) = n^^ih) + n'fih) + nl'^ih) + n\^{h) + H^f •^^^(^^^^(A;^), (67) 

-H'fih) = 70 5^ [(1 + e^(^-^^-'=^))rfl(ji, h)d^Uu h) 



Jl 

3 



^g-.2.^0-,+,)^t(^-^ + 1, k2)do{jl, k2) + h.c. , (68) 




Figure 2. Energy spectrum in tlie region around tire n — Landau level in the 
bilayer graphene with zigzag edges and a periodic boundaries orthogonal to those 
(along 62). Here we have a relatively small magnetic field (p = 1/500, with 70 = 1.0 
and 71 = 0.2. The attached numbers denote the Hall conductivity when Ep is in each 
interval between the Landau levels, which coincide with the number of edge modes 
(curves traversing different Landau levels). 



nl''{k2) = H'f^k^) : d ^ u,j, ^ Jl + 2/6, (69) 
^^(^2) = nl''ik2) -.d^t^n^ Jl + 4/6, (70) 
^r(^2) = $^[7i4(ji,M^o(ji,M 



+^3(^e^(-^2+2-'.0i+i))4(^-^, ^2) + (e-^^ + e-^2-^(^'^+t))rfl(ji + 1, fc2))wo(ji, ^^2) + h.c.Jri) 
nlf'^''''ik2) = n',^ik2) :d-^u,u^t,j,^ Jl + 2/6, (72) 
^i?'^^^(A:2) = n',^{k2) ■.d^t,u^u. (73) 

In FiglJl we have shown the spectrum of the bilayer graphene on a cyhnder. Here 
the translational symmetry along 62 makes the wave number ^2 a good quantum number, 
so that the dispersion against /c2 is depicted. The "bulk-edge correspondence", which 
dictates bulk topological features such as the quantized Hall conductivity has a one-to- 
one correspondence to edge properties, allows us to obtain the Hall conductivity of the 
system by counting the number of edge modes. 
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4.2. Zero magnetic field 

For zero magnetic field, we can make a Fourier transform along the ei as 

-^e^'^'^^Ua, da, ta{jl, k2), 

we have a fully Fourier-transformed Hamiltonian. For the bilayer we have 

d'^k 



11 



^bilayer 



n{k) = 70(1 + e-""' + e-*)(4do + uiuo) 

+liuido + lz{e-'^' + e-'^^ + e-'('=i+'=2))dUo + h.c. 



(4, 4, 4)^1 



O D 
O 
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Uo 

do 
D = 
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where 



A(/c) = 70(1 + e" 



-ik-i 



-ik'i 



73 



73 + e-'^^ + e-'''^e-'^-'). 



For the trilayer, we have 

d^k 



■^trilayer 



^Hik), 



n{k) 



7o(l + e-'^' + e-'^'){d\do + u\uo + t\to) + 



(for ABC) + -ii{t\uo + u\do) + 73(^^0 + d\uo) + h.c. 
(for ABA) + 71(^1^0 + u\do) + 73(4«o + dlu^) + h.c. 



H 



^{tluldltluldl)H, 



Q jjABC,ABA 
jjABC,ABA^ Q 



u, 
d, 

to 

Uo 

. do . 



jjABC _ 



ABA 



A 7l 

73 A 71 , L» 
73 A 

Further, we can extend the above argument to consider the ABC stacked graphene 
with a gegeral number [p) of layers. In the present representation the Hamiltonian 



A 7^ 

7i A 7i 



7^ A 



(74) 

(75) 
(76) 
(77) 

(78) 

(79) 
(80) 

(81) 

(82) 
(83) 
(84) 

(85) 



(86) 
(87) 
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simplifies into 

' O Dr. 



Dl O 



, Dp = Al + 7i Jp + 73 Jj 



pi 



where 1 is a p- dimensional unit matrix and Jp 



1 

■•• 



4-2.1. Low-energy Hamiltonians around K and K' The Dirac points (K and K' points) 
for the monolayer graphene are specified by e"*'^^ = u and e"*'^^ = o;^, = 1, w 7^ 1, at 
which A vanishes along with 73 — > 0. Around the Dirac points we can expand them as 

^ = -"^701^^, 73 = "^73^^^, ^ = 6ki + Lo6k2, (89) 

As for the multilayer graphene, gapless momenta are determined by det D's. For the 
bilayer, it is given as 



detD^^ = A^ 



7i73 



2 2 



This vanishes at four points, ^ = Q,^q,uj^q,u'^^q, where ^0 



(90) 

(7i73/7o)e'"/'. Namely, 



each Dirac point proliferates into four. 
For the trilayer, it is 

det D^^^ = det D^^^ = A^ - 2 A717;. (91) 

This gives Dirac cones at ^ = a/2^oj ^io^i ^/^io^"^ and double zero at ^ = 0. One may 
discuss the case of p— layers as well. 

The Dirac points in the chiral-symmetric system at = are specified by 
detD = 0, for which the low-energy dispersion is given by e{k) oc ib| detDl. In Fig|3l 
we have shown the low-energy dispersion along with | det D \ and arg det D for a finite 
value of 73 in the bilayer system. From this plot we can see that the chiralities of the 
four Dirac points are +1, -1,-1,-1 or -1,1,1,1. 

Now let us discuss the low-energy Hamiltonian in the absence of the trigonal 
warping (73 = 0). If we denote the number of layers by p in the ABC-stacked p-layered 
graphene (including bilayer). Dp is triangular 

z{k) 1 
z{k) 1 ■•• 



Dp{k) = 71 



det Dp{k) = A^ 



z{k) = A{k)h, 



(92) 



(93) 



where taking a continuous limit around K and K' points implies we are assuming |z| <^ 1. 
Precisely at the K and K' points, we have D^^Dp = 7idiag (0, 1, 1, ■ ■ ■); there exists 
only one low-energy mode for D^Dp while the others have 7^. Then the chiral symmetry 
(i.e.. Hp is composed of Dp as the off-diagonal block) implies that Hp has two low-energy 
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Figure 3. For the bilayer graphene in the presence of the trigonal warping (73 ^ 0), 
we show Arg det-D (left panel), contours for | det-D| (center), and detD((x the band 
dispersion; right). Here we take 70 ~ 3.1, 71 = 0.4, 73 = 0.3. 

modes, ±e(A;), with e{k) = \ det Dp{k)\/'-f^~^ = 7i|z(fc)p. Here we can take the chiral 
basis introduced in previous sections to expand the low-energy doublet as 

^. 



ipo 



DpDlij. = eV., DlDpiJo = e^iJo (94) 



ipt^o are asymptotically normalized and {\z\ ^ 1) given as 



- (_^*)P-1 - 




1 






—z 




, ^0 = 




-z* 




1 







(95) 



withijlij, = V'JV^o = {l-\z\P){l-\z\'^y = 1 + C(|2p). They are consistent with Eq.dS 



up to the errors arising from the normalization of iplip, and iplipo- Then projecting out 
the high energy sectors, we have an effective Hamiltonian for the low-energy doublet 
formed by the chiral basis ip (including the monolayer case) as 



-71 



-71) 



-(p-i) 



AP 
A*P 



(96) 



{-z)P 
{-z*)P 

As for the ABA trilayer, (D^^^y (D^^^) has two low-energy modes around the 

A 
3 



K and K' points, which implies the low-energy modes of H^^^ is spanned by a 4 





'a 
































4'' 





dimensional orthonormalized chiral basis as 

^ = (^+, ^- 

7/>V = diag(l,l,l,l) + 0(|z|), 

D,Dl^|Ji'^ = eM'\ DID,^^^ = eM'\ ^1 
where we have, up to the leading order, 

-1 



|A|, e2 = 7r'|A| 



(97) 

(98) 
(99) 



4^^^ = ^ 
V2 



= ^ 
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(1) 



1 

71 



(2) 



1 

2 



Now we have ip"^ H^^^ip 



O ^IjIDiP, 
iplD^i^, o 



end up with a simple decomposition, 



and iplDipo 



71 



z 




(100) 



we 



^ABA:eff ^ 



z* 
z 



71 





{-zf 



—z 




*\2 



Ef © iff /y2.(101) 



4-2.2. Landau levels Around the Dirac cones, we can expand the effective Hamihonian, 
as in the monolayer 



Hi ^ HI 
where Aq, 



cff 



(a ■ X)6k^ + {a ■ Y)6ky, X 



ReA^. 




ReAy 


— ImA^: 


, Y = 


—ImAy 











(102) 



§t\k=K,K'- Then we have 



eff\2 



(hcSkY, Sk = [5kx, 5ky]'E. 



6k^ 
6k„ 



X X Y 



X-X Y-X 
X-Y Y-Y 



(103) 



where = \X x Y\/fi^ is an effective "light velocity" and 5k is the averaged momentum 
near the Dirac cones with detS = (|Xp|y|2 - \X ■ Y\'^)/\X x Y\^ = 1. Note that 



(104) 



where x = sgn Im A^; A* is the chirality of the Dirac cones. 

In a magnetic field, the Hamiltonian in the continuum limit is obtained by 













X xY = 












ImA.A; 







replacing h6ka — tTq 



Pc 



cAa = —ihda — eAa where [kx, iTy 



ihcB 



with is = \/h/eB {eB > 0). As for the ABC stacked p-layered graphene (including 
monolayer and bilayer), the effective continuum Hamiltonian is obtained as the extention 
of the McCann-Fal'ko [9] as 

{Wy 

W 



h;- ^ H- ^ 



-7iJ 



-(p-1) 



(105) 



A ^ = h'^{Ax7Cx + A, 



y'^y 



(106) 



This is another derivation of the p-th layer effective Hamiltonian discussed by Koshino- 
McCann|10] for the chiral symmetric case. Here the commutation relation of H is given 

as 

[H, Ht] = 2h-Hm A,, a; = 2{ch/inyx. (107) 
Then one can define a canonical boson operator a, ([a, = 1) as 

a{ch/eB)V2 (X = +1) 



H 



{ch/eB)V2 ix 



(108) 
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Figure 4. Low-energy Landau levels in the bilayer graphene for the magnetic field 
(j) — 1/3000 as a function of 71 (with 70 = 1.0 and 73 = here). 



For the positive chirality x = +1? one has 



C\2 
P ' 







Y[Uin + ^) 



where = 71(71 ^cy2ehBY, n = a^a, and 

(a'l")^a^ = a^na = {n - l)n, {a^)^a^ = (n - 2)(n - l)n, 

a^a^y = {n + 2){n + !),■■■. 
Then we have a serise of non-equally spaced Landau levels as 

En = ±Np^/n{n - l)(n - 2) ■ ■ ■ (n - p + 1). 
The X = — 1 case follows trivially. 



(109) 

(110) 
(111) 

(112) 



4-3. Lifshitz transition in magnetic fields 

In this section, let us focus on the Lifshitz transition caused by the trigonal warping term 
73, by calculating the energy levels on a honeycomb lattice. Using the chiral symmetry, 
we can diagonalize D'^D within the low-energy sector, which suffices to discuss the 
Landau levels near E = 0. Since Z) is a sparse matrix, we can handle quite small flux (p 
numerically using the lattice Hamiltonian. 

Let us first show the energy levels of the bilayer graphene as a function of 71 in 
FigHl One can observe that larger 71 produces larger bonding- ant ibonding splitting, so 
that the low-energy sector is projected out from the high energy one. 

Figure [5] , on the other hand, displyas the Landau levels as a function of the 73. 
One can see that the Landau levels at 73 = are adiabatically connected to two classes 
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Figure 5. The Lifshitz transition caused by 73 is observed, (left) Low-energy Landau 
levels in the bilayer graphene in a magnetic field (j) = 1/3000 (70 = 1.0 and 71 = 0.5). 
The red points at 73 = denote a sequence Const. \Jn{n — 1), {n = 0, 1,2,---). 
We have also displayed the sequences oc ^/n for n = 0,1,2,3 (blue dots) and for 
n = 1,2,3 (orange), (right) The same plot for realisic 71/70 = 0.2 and much smaller 
flux = 1/30000. 

of Landau levels at 73 7^ 0, which correspond to those of two kinds of 4 Dirac cones 
with differente fermi velocity. To observe the phenomena at the moderate value of the 
flux, we take artificially large value for 71/70 = 0.5 in the left pannel of the FigsJS] In 
the right pannel, the same plot for a much smaller, realistic 71 and the magnetic field. 

5. Optical responses in bilayer and trilayer graphene 

Let us move on to discuss optical responses in bilayer and trilayer graphene in the 
quantum Hall regime. We focus on the optical longitudinal conductivity (T^s(a;), which 
descrones the optical absorption, and the optical Hall conductivity axyifjj), which 
is directly related to the Faraday rotation, since the Faraday rotation angle Qh is 
proportional to (Jxy{oj) in the quantum Hall regime as 

where Cq is the speed of light and noins) is the refractive index of the air (substrate) [TT]. 

The Faraday rotation for monolayer graphene is studied theoretically [121 [13] , and 
experimentally starts to be measured p^. Further, there are growing interests in the 
optical properties of bilayer and trilayer graphene, given that their electronic structures 
are distinct from that of monolayer graphene. As we have seen in Sec. 4, bilayer 
graphene has two parabolic bands touching at the K and K' points [T6] . Trilayer 
graphene comes with two different types of stacking, ABA and ABC, where the former 
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interband 
resonances 




Figure 6. For the bilayer graphene QHE system without the trigonal warping effect 
(73 = 0), we show (a) optical longitudinal conductivity axxi^F, w), and (b) optical Hall 
conductivity axyiep,'-^) plotted against the Fermi energy ep and the frequency lj. 

is described effectively as a mixture of monolayer and bilayer bands, while the latter 
shows a cubic band as Eq JlOSI with p = 3. Let us discuss optical responses for multilayer 
graphenes with these different lattice structures. Especially, we argue the effect of 
the presence or absence of the chiral symmetry gives a significantly different optical 
responses of n = LL in ABA trilayer graphene. 

The optical longitudinal (^^^^^(ci;)) and Hall (^^^^(ci;)) conductivities are evaluated 
from the Kubo formula as 

= — > J„ -, (113) 

^ eb- ea eb - ea - hhj - ir] 

where f{e) is the Fermi distribution, eigenenergies, rj a small energy cutoff for 
a stability of the calculation, and j'^ is the matrix element of the current operator 
j = dH/dA. 

5.1. Bilayer graphene 

We start with calculating the optical Hall conductivity using Kubo formula Eqn J113l for 
the bilayer system without the trigonal warping effect (Eq JlOSI with p = 2). The result 
for the optical longitudinal and Hall conductivities are shown in FigE]^a,b). If we label 
LLs with the Landau index n and an electron/hole band index s = ±, an intra-band 
transition (n, s) — )■ (n + 1, s) occurs around 

Wintra ~ '^bilayer = '^v'^eB/'Ji, 

with the Fermi velocity of the Dirac cone f (= c in Sec. 4), since LLs are almost equally 
spaced, unlike in the monolayer case where LL energy oc ^/n is not equally separated. 




Figure 7. For the ABA-stacked trilayer graphene QHE system (with = 5) we 

show (a) optical longitudinal conductivity axx{£F,^^), and (b) optical Hall conductivity 
f^xyi^Fi^) plotted against the Fermi energy ep and frequency uj. 



On top of this, there are inter-band transitions across the band-touching point, obeying 
a selection rule (n, —s) — )■ (n + l,s). Thus the inter-band transition occurs around 



inter 



2|eF| 



for large enough n. 

5.2. ABA trilayer graphene 

For ABA stacked trilayer, the effective Hamiltonian around Kj^/K- points is given by 
a 6 X 6 matrix (the dimension being 2 sublattices x 3 layers) as [T71 [HI [HI [20] 



Haba 








i;7r^ 







72/2 





\ 




VK 


A' 


71 








75/2 









7i 


A' 







7i 











VK 














72/2 

















v 





75/2 


7i 





VTl 


A' 


/ 



with 



TT 



ZTT. 



(114) 



a velocity V3 associated with 73 as V3 = vys/lo, ^' the on-site 



energy difference between the atoms with and without vertical bond 71, and 72(75) the 
next-nearest interlayer hoppings between Al and ^43 (51 and B3). 

First, we only retain 70,71 terms in Eq J114[ in which case the Hamiltonian 
is chiral-symmetric. In FigJ7l^a,b), the optical longitudinal and Hall conductivities 
axx{^F,^),o'xy{eF,uj) are plotted against the Fermi energy ep and frequency u. The 
magnetic field is chosen so that an interlayer hopping energy is 'ji/hwc = 5, with 
the monolayer cyclotron energy Uc 

= ^v/i = v^2eB/h. We can see a monolayer 
contribution (Dirac cyclotron frequency = Uc) and a bilayer contribution with a 
cyclotron energy, huuiayer/ V^, both of which show intra-band and inter-band transitions. 
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Figure 8. (a) Longitudinal axx{e-F,^jj) and (b)Hall axy{eF,i^) plotted against the 
Fermi energy ep and the frequency uj for a magnetic field B = IT. (c)A diagram 
indicating allowed resonances in axy- 



For moderate magnetic fields ~ 1 T, cyclotron frequency for monolayer is much larger 
than that for bilayer Uc ^ c^biiayer- In this case, a^x is an even function with respect to 
ep = 0, while axy is odd due to the electron-hole symmetry, a consequence of the chiral 
symmetry. The jump of axy at e^? = is related to the chiral zero modes. 

Now, if we include all the hopping terms in Eq J114l the chiral symmetry is broken, 
and we no longer have the chiral protection for zero modes (three zero energy LLs). 
We have then a massive Dirac band plus gapped bilayer bands with energy shifts. As 
for hopping parameters, we adopt the values for graphite, 70 = 3.2 eV, 71 = 0.39 eV, 
73 = 0.32 eV, 72=-0.020eV, 75=0.038eV, A'=0.050eV [211 [7]. Figures [7(a,b) depict the 
result for the optical longitudinal axxi^F,^^) and Hall conductivity axy{eF,uj) plotted 
against the Fermi energy ep and frequency u for ABA stacked trilayer graphene QHE 
system. We see contributions from monolayer-like Dirac LLs (labeled with M) and those 
from bilayer LLs (B), both of which exhibit intra-band and inter-band transitions. Since 
Dirac cone is massive due to the chiral symmetry breaking and the zero-energy LL for 
the monolayer band (MO) is situated at the bottom of the conduction band for Kj^ valley 
and the top of the valence band for K_, MO — M(l, +) resonance occurs at an energy 
lower than M(l, — ) MO for and vice versa for K^. A cancellation of resonances 
in axy, due to the opposite signs in current matrices, occurs between M(l, — ) — )• MO 
for and MO — )■ M(l, +) for for a region of Fermi energy between MO{Kj^) and 
MO{K_.), while this is not the case with axx- For bilayer contributions, satellites appear 
as B{n, ±) Bin + 1 + 3m, ±) and B[n, ±) ^ B{n + 2 + 3m, ±), since the trigonal 
warping term mixes n LLs and n + 3 LLs [22] . The resonance frequency for intra-band 
transition within the conduction band is larger than those within the valence band, 
which is a consequence of the electron-hole asymmetry in the bilayer bands, triggered 
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inter-band 




Figure 9. For ABC-stacked trilayer graphene QHE system (with ji/hujc = 5) we 
show (a) optical longitudinal conductivity (^^.^.(ei?, cj), and (b) optical Hall conductivity 
(Jxyi^FjUj) plotted against the Fermi energy ep and frequency uj. 

by a breaking of the chiral symmetry. A deviation in the cyclotron mass for electron 
and hole bands prevents a complete cancellation between B{n, — ) — t- B{n + 1, +) and 
B{n + 1,—) — 7- B{n,+) transitions, which results in small interband transitions in a 
wide region of Fermi energy. 

5.3. ABC trilayer graphene 

The low energy effective Hamiltonian of ABC stacked trilayer graphene in a 2 by 2 
matrix Eq JlOSl is a cubic form in the momentum, if we neglect hopping terms but 7o, 7i- 
For this effective Hamiltonian, LL energy (Eq J112p shows a magnetic-field dependence 

3 1 

oc 5 2, resulting in a smaller LL spacing compared to the single-layer LL oc B2 and 
bilayer LL oc B for weak magnetic fields. 

Now we turn to a result for the optical conductivities in ABC trilayer for an 
interlayer hopping energy 'fi/huc = 5, with monolayer cyclotron frequency Uc = 
v\j2eB/h. Figures |9]|^a,b) show the optical longitudinal and Hall conductivity axyiep, 00) 
plotted against the Fermi energy ep and frequency u for ABC stacked trilayer graphene 
QHE system. In contrast to ABA stacking (Fig I?]), we only see a single sequence of 
intra-band and inter-band transitions with a much smaller cyclotron energy than 

3 

due to the dependence on magnetic fields oc B2 . The inter-band transition occurs at 
hjjj 2eF for the same reason as in the bilayer, while the intra-band transition energy 
grows with increasing LL index n as 

which explains why the intra-band resonance energy increases with n in FiglHI^d), 
while it decreases with n in the case of the monolayer graphene [13]. We can then 
predict in general that, for ABC p-layered graphene, the intra-band transition occurs 
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at and the inter-band transition at u ~ 2ep. Thus we end up with the 

intra-band transition energies that exhibit different behaviors with Landau index n for 
monolayer (oc n~^/^), bilayer (constant) and ABC trilayer (oc n^/^) graphenes, while the 
inter-band transition energies are qualitatively the same with ~ 2ei?. 
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